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Abstract 

Small polaron behavior in a two dimensional honeycomb net is studied by 
applying the strong coupling perturbative method to the Holstein molecular 
crystal model. We find that small optical polarons can be mobile also if the 
electrons are strongly coupled to the lattice. Before the polarons localize and 
become very heavy, there is infact a window of e-ph couplings in which the 
polarons are small and have masses of order ~ 5 — 50 times the bare band mass 
according to the value of the adiabaticity parameter. The 2D honeycomb net 
favors the mobility of small optical polarons in comparison with the square 
lattice. 
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I.INTRODUCTION 



The transport properties of real systems are strongly affected by the presence of non lin- 
ear potentials in the lattice Nonlinearities may arise because of embedded impurities 
in the host lattice which favor trapping of the charge carriers [ffl-ffl] in any dimensionality || 
or, nonlinearities may be intrinsic to the system and driven by the electron-lattice interac- 
tion as assumed in the Molecular Crystal Model first proposed by Holstein [[J in the form 
of a discrete nonlinear Schrodinger equation for electrons coupled to harmonic phonons. 
Several theoretical methods have been developed in the last decades [|TT|-|iT)[ to analyse the 
ground state and finite temperature properties of the unit comprising the electron plus the 
surrounding local lattice deformation, namely the polaron. Central to these investigations 
is the concept of self-trapping traditionally denoting a transition between an infinite radius 
state at weak electron-phonon coupling and a finite radius polaron at strong e-ph coupling. 
The narrowing of the polaron band and the abrupt increase of the polaron effective mass ver- 
sus e-ph coupling are the classical and related indicators of the transition event which may 
occur or not according to the degree of adiabaticity and the dimensionality of the system. 
When the characteristic phonon energy u becomes larger than the electronic bandwidth 
the antiadiabatic regime is attained. In this case it is generally accepted that the polaron 
wave function spreads over a few lattice sites (small polaron) with the polaron mass being a 
smooth function of the e-ph coupling. Instead, the polaron self-traps in the adiabatic regime 
and there is growing evidence that this event takes place in any dimensionality [|lf|,|ll].[2(3 . 

When the lattice polarizations is confined to one or a few unit cells the carrier feels the 
details of the local structure: recent generalizations of the Holstein model have shown that 



the inclusion of on site lattice anharmonicity can substantially modify the size 21 and the 



mass 



22] of the polaronic quasiparticle. While these findings could contribute to locate with 



more accuracy the self-trapping event in parameter space it is still unclear whether and how 
the transition depends on the lattice structure. To address this problem we focus here on the 
polaron mass renormalization in a two dimensional honeycomb net which can be viewed as 
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a triangular Bravais lattice with a bases of diatomic molecules at the vertexes. The Holstein 
model is briefly reviewed and the results are discussed in Section II. The conclusions are 
drawn in Section III. 



II. THE HOLSTEIN MODEL HAMILTONIAN 

The Hamiltonian for the single electron in the Holstein model reads 

h = -tj24 c j + 9j2 c l c i( a i + a l) + J2 uj k a l a k (Q- 1 ) 

i^j i k 

where the dimension dependence explicitly appears in the momentum space Hamiltonian 
for the harmonic lattice vibrations, cj ( q ) creates (destroys) a tight binding electron at the 
i site and t is the first neighbor hopping integral related to the bare electron half bandwidth 
D by D = zt, z being the coordination number, ( ) creates (destroys) a k-phonon 
with frequency u^. g is the overall electron-phonon coupling constant. 

In the strong coupling regime the Lang-Firsov approach [^3| is reliable |2^ , |25|| and the 



polaron mass m* can be obtained via perturbative method. In d dimensions the ratio 
between m* and the bare band mass m is p0|: 



fm*\ exp(g 2 ) 

^m$' d 1 + z 2 texp(—g 2 ) f '(g 2 ) / \huj) 

-2 ^9 sr^ ■ 2 sr^ \-2 
9 = jj- 2^ sin — }^ (hub) 

The series expansion in the last of eqs. (2) reflects the fact that the second order polaron 
self energy comprises the emission and absorption of an arbitrary number of phonons hence 
it is a sum over an infinite number of diagrams each having s phonons between the two 
interaction vertexes. The second order of perturbative theory also introduces the effect 
of the adiabaticity parameter zt/hu on m* which, in general, depends on dimensionality 
through: i) g 2 oc d, ii) the first neighbors number z, iii) the Brillouin zone sums and iv) the 
features of the phonon spectrum. 



We take a 2D honeycomb net equivalent to a triangular lattice with a two points basis. 
Each lattice site is a diatomic molecule with coordination number z — 6. Hence, the 
phonon spectrum has both acoustical and optical branches whose analytical expressions can 
be deduced by a force constants parametrization scheme: 

/3 + 3 7 



^(k Xl ky) = i — • yJP* + J*H(k x , ky) + /3jG{k x , ky) 

Hyk X , ky) = 3 + l{C X Cy + S X Sy + C^ X Cy + S^ X Sy + c 2 

G{k x ,k y ) = 2(2c x c y + c 2x ) (0.3) 

with: c x = cos(k x \/3a/2), c y = cos(k y 3a/2), C2 X = cos{k x \/3a) ) c 3x = cos(k x 3y/3a/2), 
s x = sin(k x \/3a/2), s y = sin(k y 3a/2), s 3x = sin(k x 3\/3a/2). a is the lattice constant and 
M is the reduced molecular mass. (3 and 7 are the intra- and mier-molecular force constants 
respectively in terms of which one defines: Uq = 2(3 /M, io\ = 7/M and the zone center 



optical frequency uj = \JujI + zuj\. 

Previous investigations of the Holstein Hamiltonian |2(| have shown that the intermolec- 
ular forces have to be sufficiently strong in order to predict the correct polaron bandwidth 
trend versus dimensionality. Thus the ground state properties of the Holstein model essen- 
tially depend on the strength of u)\ which should be of order ~ cuo/2- When this condition 
is fulfilled the polaron mass turns out to be substantially dimension independent. Larger uoi 
values are admitted in the model although they encounter the obvious upper bound uj\ < uj 
in a molecular lattice. This result (which has been proven in a large portion of parameter 
space ranging from fully adiabatic to antiadiabatic conditions) introduces a novel feature 
in the polaron landscape corroborated by Monte Carlo simulations | 27fl , density matrix 
renormalization-group studies and variational approaches. 



Two quantities play a central role in polaron theory. The first, defined by 

Ng 2 

iV being the number of molecular sites, represents the ratio PS[ between the polaron binding 



energy and the electronic half bandwidth. It yields the energetical gain due to small polaron 



formation with respect to the bare electronic state. The second, defined by 

Ng 

a = — 

Ek n^k 

, measures the lattice deformation associated with the quasiparticle formation. While in 
adiabatic systems the condition A > 1 signals the existence of the small polaron state, in 



antiadiabatic systems a > 1 is a more restrictive condition for small polaron formation |29 

Recent analysis |3O|-0^] on the mobility of small polarons also in conjunction with models 
on polaronic high T c superconductivity [[33H 3511 have led to reconsider the concept of self- 
trapped state which, although being intimately related to the small size of the quasiparticle, 
is not synonymous with small polaron state. Thus, if the narrowing of the polaron bandwidth 
(induced by an increasing e-ph coupling) marks the onset of the transition between large 
and small states still there is a range of g values for which the polaron, although spread over 
a few lattice sites only, is not trapped and retains mobility properties. The self-trapping 
event can be instead associated with a rapid but continuos effective mass increase which is 
precisely located by looking at the curvature of the logarithm of the effective mass versus g. 
In our view this method, beyond embodying the full physical significance of the transition 
process, offers a simple criterion to select a "critical g value" as an inflection point either 
in the logarithm of the effective mass or in its first derivative. While the former case would 
identify the point of most rapid increase of the effective mass, the occurence of the latter case 
distinguishes a peculiar point in the mass versus g plot although the concavity-convexity 
change is absent. 

In Fig. 1(a), four plots of the polaron mass (in units of the bare band mass) as a function 
of the g coupling are reported on while the corresponding curvatures of the logarithm of the 
mass are shown in Fig. 1(b). High optical phonon frequencies are assumed. Our selected plots 
range from an extreme adiabatic (t = 200meV) to a moderately antiadiabatic (t = lOmeV) 
regime. At a fixed g the antiadiabatic polaron is always heavier than the adiabatic polarons 
but the peculiar points, resolved as the minima of the curves displayed in Fig. 1(b), lie at 
decreasing g values by decreasing the degree of adiabaticity. Then, an antiadiabatic polaron 



may self-trap already at g ~ 2.7 weighing m e // ~ 14 while an extreme adiabatic polaron 
self-traps only at g ~ 3.8 with wi e // — 123. The correctness of our perturbative approach 
is monitored by the A and a values which are larger than one in all the displayed points 
consistently with the assumption on the existence of small polarons. In the intermediate 
adiabatic cases t = lOOmeV and t = 50meV we find m e ff ~ 73 at g ~ 3.5 and m e ff ~ 40 at 
g ~ 3.3 respectively. Thus, a 2D honeycomb net seems to sustain mobile adiabatic polarons 
in a window of strong coupling regimes approximately defined by 2.5 < g < 3.3. In Figs. 2, 
the phonon frequencies are still high, although much reduced with respect to Figs.l. As a 
main effect the polaron masses are roughly doubled while the peculiar minima of the second 
derivative of the logarithm of the effective mass (Fig.2(b)) do not shift substantially versus 
g with respect to the corresponding cases in Fig. 1(b). To emphasize the role of the lattice 
structure we have compared the adiabatic polaron (with t = lOOmeV) in the present 2D 
honeycomb net with the previously investigated |2(| square lattice. As an example, provided 
that: i) the same overall e-ph coupling g = 3.38 is taken, ii) the same values of intra- and 
inter-molecular force constants are assumed (uq = 50meV and u>i = 25meV) and these 
values are consistent with the strong coupling perturbative method, we find m e ff ~ 140 for 
the honeycomb net polaron against m e ff ~ 1200 for the square lattice. Also the intermediate 
adiabatic polaron (t = 50meV) behaves in a similar manner being m e ff ~ 200 in the 
honeycomb net polaron against m e ff ~ 1700 in the square lattice assuming the same input 
parameters as above. We have also considered the effect of the acoustical branch of the 
phonon spectrum on the polaron properties. In all cases, with different degree of adiabaticity, 
the minima of the second derivative (Fig.3(b)) occur at much lower g values than for optical 
polarons. However, no physical meaning can be attached to these values since they are 
well outside the range of validity of the Lang-Firsov based perturbative method. When the 
method holds (g > 2.5) no peculiar point can be resolved in the polaron behavior versus g 
which is anyway characterized by a huge mass renormalization. Figs. (3) have been reported 
also to point out that the occurence of distinctive features in the mass or mass derivative 
curves, far from being a cogent criterion for the self-trapping event, just indicates a trend 
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which needs to be corroborated by the analysis of other independent quantities. 
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FIGURES 



FIG. 1. (a) Polaron masses (in units of bare band electron mass) versus g in adiabatic and 
antiadiabatic regimes with high frequency optical phonons: u = lhSmeV . (b) Second derivative 
of the logarithm of the effective mass with respect to g. 
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FIG. 2. (a) Polaron masses (in units of bare band electron mass) versus g in adiabatic and 
antiadiabatic regimes. The characteristic frequency of the optical phonons is: u = 79meV. (b) 
Second derivative of the logarithm of the effective mass with respect to g. 
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FIG. 3. (a) Polaron masses (in units of bare band electron mass) versus g in three adiabatic 
regimes. The electron couples here to the acoustical phonon branch, (b) Second derivative of the 
logarithm of the effective mass with respect to g. 

IV. CONCLUSIONS 

We have developed a perturbative study of the molecular crystal model assuming the 
existence of strong eiectron-phonon coupling conditions which favor the formation of small 
polarons. The Lang-Firsov method permits to calculate the mass renormalization for spe- 
cific structures once the phonon spectrum is known. Rather than applying the model to real 
systems as previously done for simple lattices, we have examined whether alternative struc- 
tures as the 2D honeycomb net may host polarons which are both small and mobile. Infact 
we have found that adiabatic small polarons are lighter by roughly a factor eight than in the 
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square lattice once the same input parameters are assumed. Adiabatic polarons have been 
studied as a function of the e-ph coupling and distinctive points in the mass behavior (versus 
g) have been selected and put in relation with a possible occurence of the self-trapping event. 
Although not ultimative this criterion seems plausible provided that such "self-trapping g 
points" are obviously well within the range of applicability of the Lang-Firsov scheme. We 
find that small optical polarons in the honeycomb net can self-trap if the electron-lattice 
system couples in the range 3 < g < 4. The exact location of the transition depends on 
the adiabaticity parameter with intermediate adiabatic polarons lying in the lower portion 
of that range. Our results suggest that adiabatic optical polarons can be mobile in the 
honeycomb net although a strong coupling regime holds. The small polaron effective mass 
is of order ~ 5 — 50 times the bare band mass before the self-trapping point is attained, with 
more adiabatic polarons being lighter once the e-ph coupling is fixed. 
This work has been done in Trequanda (Siena). 
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